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Abstract 

We consider the correlation between the binding energies of the triton and the a-particle which 
is empirically observed in calculations employing different phenomenological nucleon-nucleon in- 
teractions. Using an effective quantum mechanics approach for short-range interactions with large 
scattering length \a\ S> £, where £ is the natural low-energy length scale, we construct the effective 
interaction potential at leading order in £/\a\. In order to renormalize the four-nucleon system, it 
is sufficient to include a 5'f7(4)-symmetric one-parameter three-nucleon interaction in addition to 
the S-wave nucleon-nucleon interactions. The absence of a four-nucleon force at this order explains 
the empirically observed correlation between the binding energies of the triton and the a-particle. 
We calculate this correlation and obtain a prediction for the a-particle binding energy. Corrections 
to our results are suppressed by £/\a\. 
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The scattering of two particles with short-range interactions at sufficiently low energy 
is determined by their S-wave scattering length a. Generically, the scattering length a is 
comparable to the typical low-energy length scale of the system £ which, for short-range 
interactions, is of the order of the range of the potential. In special cases, however, the 
scattering length can be much larger, \a\ 3> £, due to a fine-tuning of the parameters in 
the underlying interaction potential. Such a fine-tuning can be accidental as in the case of 
two nucleons at low energy, or it can be controlled experimentally by varying an external 
parameter such as a magnetic field in the case of Feshbach resonances for atomic systems. 

Few-body systems with large scattering length are particularly interesting because such 
systems show many universal properties. The simplest example for positive a is the existence 
of a shallow bound state (dimer) with binding energy 



where m is the mass of the particles. The separation of scales between \a\ and £ also 
manifests itself in few-body systems with more than two particles. The most prominent 
example is the Efimov effect: the occurence of a geometric spectrum of three-body bound 
states in the limit a — > ±00 1]. In this limit, the ratio of the binding energies of successive 
Efimov states approaches the universal number 515.03.. . In order to fully specify the bound 
state spectrum, however, the absolute energy of one Efimov state has to be specified and 
this introduces one three-body parameter, L 3 , which will be discussed in more detail below. 
As a consequence, a three-body system with a large but finite scattering length can be 
characterized by specifying the values of a and L 3 . 

There has been a considerable interest in physical systems with large scattering length 
recently. This interest was stimulated by the experimental realization of Feshbach resonances 
with trapped atoms and by the effective field theory program for nuclear physics. (See, e.g., 
Refs. 0, H[ and references therein.) The best known example of a nuclear system with a 
large scattering length is the two-nucleon (NN) system. There are two independent S-wave 
scattering lengths that govern the low-energy scattering of nucleons. The scattering lengths 
a s = —23.5 fm and a t = 5.42 fm describe NN scattering in the spin-singlet ( 1 5'o) and spin- 
triplet ( 3 Si) channels, respectively. Both scattering lengths are significantly larger than the 
natural low-energy length scale £ ~ Tijim^c) ~ 1.4 fm, while the effective ranges are of the 
same order as £. As a consequence, an expansion in £/\a\ is useful. 

A peculiar universal feature of the three-nucleon system is the Phillips line 0. If the 
predictions of different nucleon-nucleon potentials for the triton binding energy B t and the 
spin-doublet neutron-deuteron scattering length a£j are plotted against each other, they 

fall close to a line. This correlation between B t and a^l is called the Phillips line and can 
not be understood in conventional potential models. However, it immediately follows from 
universality if the large NN scattering lengths are exploited within an expansion in £/\a\ 
If corrections of order £/\a\ are neglected, all low-energy 3-nucleon observables depend 
only on the spin-singlet and spin-triplet scattering lengths a s and a t and the three-body 
parameter L 3 . Since the NN potentials reproduce the scattering phase shifts, they all have 
the same scattering lengths. The short distance part of the potentials which is encoded in 
the three-body parameter L 3 , however, is not constrained by the phase shifts and in general 
is different for each potential. The different potential model calculations must therefore 
fall close to a line which is parametrized by the parameter L3. A similar universal feature 
of the four-nucleon system is the Tjon line: an approximately linear correlation between 




a > , 



(1) 



2 



the triton binding energy B t and the binding energy of the a-particle B a . This correlation 
was discovered by Tjon p] using simple separable interactions, but also holds for modern 
phenomenological potentials |8[. The origin of this correlation has not been explained. 
However, it would immediately follow if the four-nucleon system could be renormalized to 
leading order in l/\a\ without specifying a four-body parameter. The purpose of this paper 
is to study the four-nucleon system in an effective quantum mechanics approach (see, e.g., 
Ref. to leading order in i/\a\ in order to understand the origin of the Tjon line. We will 
show that no new four-body parameter is required to renormalize the four-nucleon system 
and calculate the binding energy of the a-particle to leading order in £/\a\. 

Effective theories are a convenient tool to calculate observables of low-energy systems 
to high precision in a systematic fashion. They are ideally suited to exploit a separation 
of scales such as the one between \a\ and i. At sufficiently low energy, every interaction 
can be considered pointlike. One can then use a very general effective theory to describe 
the universal low-energy properties of the system in which observables are computed with a 
controlled expansion in £/\a\. Few-body systems with large scattering length are of special 
interest: It is known that in the three-body sector with short-range interactions a three- 
body force at leading order is needed to renormalize the theory [10(. This implies that at 
leading order in l/\a\, the properties of the three-body system with large scattering length 
are not determined by two-body data alone and one piece of three-body information (such 
as the three-body binding energy) is required as well. Recently we considered the four-boson 
system with short-range interactions. We showed that at leading order no four-body force 
is needed for consistent renormalization and applied the theory to atomic 4 He clusters [ill ]. 

In this letter, we extend this work to the four-nucleon system. There has been a large 
interest recently in applying effective field theory methods to the few-nucleon problem. (See, 
e.g., Refs. 0,0, k| and references therein.) The four-nucleon problem is more complicated 



than the four-boson problem due to spin and isospin. Both the 1 5'o and 3 Si NN scattering 
lengths a s and a t contribute at leading order in i/\a\. However, as in the case of bosons, 
only one three-body parameter L3 enters for the a particle. The parameter L3 determines 
the S , f/(4]-symmetric three-nucleon force that is required to renormalize the three-nucleon 
problem p]. In the following, we will show that no four-body parameter is required to 
renormalize the nuclear four-nucleon system and that the binding energy of the a-particle 
can be described to 10% accuracy at leading order in i/\a\. Furthermore, we will calculate 
the Tjon line at leading order in £/\a\. 

The effective low-energy interaction potential generated by short-range contact interac- 
tions can be written down in a momentum expansion. In the S-wave sector of the two-nucleon 
system, it takes the general form 

(k'|V|k) = V s X s 2 <?(k')<7(k) + P t A^(k')<?(k) + . . . , (2) 

where the dots indicate momentum dependent terms that are higher order in £/\a\ and 
V s and Vt project onto the 1 S'o and 3 Si partial waves, respectively. Because of Galilean 
invariance, the interaction can only depend on the relative momenta of the incoming and 
outgoing particles k and k'. The regulator function (u\g) = g(u) = exp(— u 2 /A 2 ) cuts off 
the contribution of momentum states with u ^> A. The cutoff parameter A is arbitrary 
and in the end all observables should be independent of A. The interactions in Eq. (j2J) are 
separable and thus, the two-body problem for each partial wave can be solved analytically. 
The two-body t-matrix can be written as t t>s (E) = \g)r ttS (E)(g\, where E denotes the energy. 
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The two-body propagator Tt )S (E) is given by: 



T t,s(E) = [1/AST-47T J dqq'j^}' 1 • (3) 

Here and in the following we set h = m = 1. The coupling constants A| and A2 can be fixed 
by demanding that the triplet and singlet scattering length a t and and a s are reproduced 
correctly by the corresponding t-matrices. 

The properties of the triton are determined by the Faddeev equations. To leading order 
in £/\a\, all internal orbital angular momenta can be set to zero. Thus, the triton spin 1/2 
is built up from the spins of the nucleon only. The triton wave function can be decomposed 
into Faddeev components. Since we are mainly interested in the binding energies and not in 
the wave functions, we can eliminate all but one of the components and obtain an equation 
for remaining component ip : ip = GotPip + GotGot 3 (l + P)ip , where Go is the free three- 
particle propagator and P = P13P23 + -P12-P23 is a permutation operator that generates the 
omitted Faddeev components; P^- exchanges particles i and j. The auxilliary quantity t 3 is 
the solution of a Lippmann-Schwinger equation with a S'?7(4)-symmetric three-body contact 
interaction 

V3 = Va^zKux, u 2 )h(u' 1 , u' 2 ) (4) 

only. Here V a denotes the projector on the total antisymmetric three-body state with total 
spin S = 1/2 and total isospin T = 1/2 as for example given in 0]. The regulator function 
h(ui,U2) = exp(— (u\ + I^D/A 2 ) is defined in terms of the familiar Jacobi momenta of the 
three-body system. The three-body coupling A3 must be determined from a three-body 
observable. This interaction is required in order to renormalize the three-body system and 
achieve independence of the cutoff parameter A. Its renormalization group behavior is 
governed by a limit cycle jy, |llj . For large values of A the running of the coupling constant 
is described by 

_c_ sin(g ln(A/L 3 ) - arctan(l/s )) 

31 ~~ A 4 sin(s ln(A/L 3 ) +arctan(l/s )) ' 

where c ~ 0.016 is a normalization constant and sq ~ 1.00624 is a transcendental number 
that determines the period of the limit cycle. If the cutoff A is multiplied by a factor 
exp(n7r/so) ~ (22.7) n with n an integer, the three-body coupling A3 is unchanged. The 
three-body parameter L 3 can be determined directly from observable quantities like the 
triton binding energy B t . 

Having expressed all relevant coupling constants in terms of physical observables, we are 
now ready to compute the binding energies of the four-nucleon system. 

The binding energies of a four-body system of nucleons can be determined with the 
Yakubovsky equations. The four-body wavefunction \1/ is first decomposed into Faddeev 
components, followed by a second decomposition into so-called Yakubovsky components. In 
the case of nucleons one ends up with two components if) a and ipB- In the following, we will 
consider the a-particle with total angular momentum J = and total isospin T = 0. The 
Yakubovsky equations can be written as (l5| : 

"ipA = G tP [(1-^34)^ + ^5] +^(l + G t 12 )G V 3 ^, 

4> B = G tP [(1-^34)^ + ^1 , (6) 
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where Go is the free four-particle propagator and \I/ is the full four-body wave function 
that can be reconstructed from ipA and ipB- We solve Eqs. (jHJ) numerically in momentum 
space. For details of the solution and structure of the above equations, we refer the reader to 
Refs. [IIElI. The two-body couplings A2 and A| are fixed by matching Eq. (jHJ) to reproduce 
the triplet and singlet scattering lengths a t and a s , respectively. Instead of a t one can also 
use the deuteron binding energy B d = 2.22 MeV together with Eq. (Q). The difference 
between these two methods is higher order in £/\a\ and thus gives a naive estimate of the 
error made in the order considered. The three-body parameter L3 which is implicit in the 
three-body interaction V3, is fixed from the triton binding energy B t = 8.48 MeV. We find 
L 3 = 21.2 fm" 1 (L 3 = 19.8 fm" 1 ) if a t (Bd) are used as input. 

We now calculate the a-particle binding energy B a for different values of the cutoff 
parameter A. We vary A in the region where the three-nucleon system has exactly one 
bound state. For values of A below 4 fm -1 there is a strong cutoff dependence. For values 
in the range 8 fm" 1 < A < 10 fm -1 , however, our results are numerically stable and 
independent of A to within 5%. This variation is considerably smaller than the errors from 
higher orders in £/\a\ of about 30%. The existence of the plateau for A > 8 fm -1 provides 
sufficient evidence that the renormalization of the four-nucleon system at leading order in 
£/\a\ does not require a four-nucleon force. Renormalization of the the three-nucleon system 
automatically guarantees renormalization of the four-nucleon system. 

For the a-particle binding energy, we find B a = 29.5 MeV. If the deuteron binding energy 
B d is used as input instead of the triplet scattering length a t , we obtain B a = 26.9 MeV. This 
variation is consistent with the expected 30% accuracy of a leading order calculation in £/\a\. 
Our results agree with the (Coulomb corrected) experimental value B^ p = (29.0±0.1) MeV 
to within 10%. We conclude that the ground state energy of the a-particle can be described 
by an effective theory with short-range interactions only. Of course, a more refined analysis 
should include Coulomb corrections. However, the size of these corrections is expected to be 
smaller than the NLO contribution. Furthermore, corrections due to isospin violation will 
appear naturally at higher order within the effective theory. 

The empirically observed correlation between the binding energies of the triton and the 
a-particle (the Tjon line) follows immediately from the absence of a four-nucleon force at 
leading order in £/\a\. The Tjon line is generated by variation of the three-body parameter 
L 3 , which is different for different AW potentials. In Fig. Q we show our result for the 
Tjon line with a s and B d as input (lower line) and a s and a t as two-body input (upper line). 
Both lines generate a band that gives a naive estimate of higher order corrections in £/\a\. 
We also show some calculations using phenomenological potentials ||| and a chiral EFT 
potential with explicit pions 17, 3|- The Tjon line is a general feature of few-body systems 
with large scattering length. As a consequence, all calculations with interactions that give a 
large scattering length should lie close to this line. Different short- distance physics and/or 
cutoff dependence should only move the results along the Tjon line. This can for example 
be observed in the NLO results with the chiral potential indicated by the squares in Fig. Q 
or in the few-body calculations with the low-momentum AW potential V\ ow k carried out in 
Ref. 0| . The V\ ow k potential is obtained from realistic NN interactions by intergrating 
out high-momentum modes above a cutoff A but leaving two-body observables (such as the 
large scattering lengths) unchanged. The results of few-body calculations with V\ ow k alone 
are not independent of A but lie all close to the Tjon line (cf. Fig. 2 in Ref. 0]). (For a 
different approach to calculate few-nucleon systems based on a low momentum potential, 
see Ref. |20[ and references therein.) A similar conclusion holds for the study with various 
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FIG. 1: The correlation between the binding energies of the triton and the a-particle (the Tjon 
line). The lower (upper) line shows our leading order result using a s and (a s and at) as two-body 
input. The grey dots and triangles show various calculations using phenomenological potentials 
without or including three-nucleon forces, respectively. The squares show the results of chiral 
EFT at NLO for different cutoffs while the diamond shows the N 2 LO result. The cross shows the 
experimental point. 



nonlocal potentials in |2l|. Our result for the Tjon line to leading order in £/\a\ is shifted 
from the empirical values because we neglect finite range effects. Furthermore, the empirical 
Tjon line becomes a narrow band because of higher order three-nucleon force effects. 

In summary, we have computed the a-particle binding energy at leading order in the 
effective theory with contact interactions. Our results are in good agreement with the 
experimental value and agree within the expected error bounds given by an estimate of 
higher order contributions. At leading order no four-body force is needed for renormalization 
and thus, the Tjon line can be understood as a general feature of low-energy systems with 
large scattering length. As a consequence, all calculations with interactions that give a large 
scattering length should lie close to this line. Finally, further effort should be devoted to 
the computation of scattering observables and the inclusion of effective range corrections. 
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